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CHARACTERIZATION AND DESIGN MECHANICS FOR 
FIBER-REINFORCED METALS 
by Christos C. Chamis 
Lewis Research Center 

SUMMARY 

Certain aspects of fiber composite mechanics are utilized to characterize fiber- 
reinforced metals in the elastic range of the fiber when the properties of their constit- 
uents are available and the fabrication processes are known. The concept is illustrated 
by selecting a boron aluminum unidirectional fiber composite system. This system is 
characterized from the geometry, thermoelastic response, and strength behavior view- 
points. Some comparisons with limited experimental data are made for correlation 
purposes. Theoretical results of the boron aluminum composite are presented to illus- 
trate the trends in the thermoelastic and strength properties and to point out difficulties 
in experimental work and design. The detrimental effects of voids (porosity) on some of 
these properties are also illustrated. Failure envelopes for combined stress and off- 
axis loading are constructed as an aid to designers. Some results for angle-ply com- 
posites are presented. References dealing with other aspects of composite mechanics, 
not covered herein, are cited. Other design factors such as ply thickness, density, 
delamination, and buckling parameters are discussed, and design curves for these fac- 
tors are plotted. 

The composite mechanics proposed herein provides a useful tool for the complete 
linear characterization and generation of design data of fiber-reinforced metals. Areas 
are pointed out where additional experimental data and further research are needed. 


INTRODUCTION 

The technology of fiber-reinforced metals has advanced to the point where actual 
hardware components are being fabricated from these materials. References 1 to 5 
provide a thorough review of the state-of-the-art of fiber- reinforced metals. Several 
of their underlying principles, their various fabrication processes, and some of their 
peculiarities are described in these references. Some systems that have been success- 


fully fabricated and hold promise for structural components are carbon-filament- 
reinforced cobalt and nickel (ref. 6), beryllium-fiber-reinforced aluminum (ref. 7), 
refractory-fiber- reinforced superalloys (refs. 8 and 9), boron- fiber- reinforced titanium 
(ref. 10), graphite-fiber- reinforced nickel (ref. 11), and boron-fiber- reinforced alum- 
inum (ref. 12). Studies on the compatibility of fibers with metal matrices are reported 
in reference 13. The joining of metallic matrix composites is described in reference 14. 
One specific application of composites to blading for gas turbine engines is illustrated 
in reference 15. Additional references on fabrication processes are cited in refer- 
ence 3. 

The available work to date is significant in many respects and has resulted in the ac- 
complishment of several research objectives such as determination of several possible 
composite systems, demonstration of several possible fabrication processes, and reso- 
lution of constituent compatibility, elementary mechanics, and experiments for constit- 
uent load transfer. The experimental and theoretical work which leads to efficient 
component design from fiber-reinforced metals is still in its infant stage. Herein, an 
important portion of the theoretical aspect of the information domain is treated. The 
theoretical treatment is supplemented with recommendations on experimental and addi- 
tional theoretical work which will help complete the major part of this domain. 

The theoretical portion of the information domain contains theories for structural 
and stress analysis, representative macro- and micromechanics, mathematical models, 
and efficient optimization algorithms. Structural and stress analysis and optimum de- 
sign of fiber composites, in general, require knowledge of several thermoelastic and 
strength properties and certain other design factors. Composite properties depend 
strongly on the quantities of constituent materials present and the angle between the 
loading and the fiber directions. The large number of possible constituent combinations 
and fiber orientations in a composite make it considerably difficult, if not impossible, 
to determine experimentally the thermoelastic and strength properties required in de- 
sign and analysis. The other alternative is to construct micromechanics theories 
(ref. 16) to predict several or all of these properties for unidirectional fiber compos- 
ites (UFC or plies) (fig. 1). Then lamination theory can be used for composite proper- 
ties in general. The thermoelastic and strength behavior of unidirectional composites 
is influenced by several factors, some of which are quite difficult to incorporate in 
mathematical models (ref. 16). Most of these difficulties are bypassed by developing 
semiempirical theories. 

The semiempirical theories proposed in references 17 and 18 have been success- 
fully applied in predicting all the thermoelastic and strength properties of fiber, 
nonmetallic-matrix composites which exhibit linear behavior. Herein, these theories 
are applied to predict all of the thermoelastic and strength properties of fiber, metallic- 
matrix composites in the linear fiber regime. 
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Figure 1. - Boron aluminum ply geometry. 


The parameters necessary to describe completely the thermoelastic and strength 
behavior of unidirectional fiber composites (fig. 1) consist of seven thermal properties 
(three thermal coefficients of expansion, three conductivities, and one heat capacity), 
nine elastic properties (three normal or extensional moduli, three shear moduli, and 
three Poisson ratios), five uniaxial strengths (longitudinal tensile, longitudinal compres- 
sive, transverse tensile, transverse compressive, and intralaminar shear), and also a 
combined stress strength criterion. In addition, factors such as ply thickness, inter- 
fiber spacing, density, ply weight per unit area, transverse shear delamination crite- 
rion, buckling coefficients, and weight considerations are important in practical design 
and tradeoff studies. 

The semiempirical theories described in references 17 and 18 are applied to some 
available experimental results (refs. 19 and 20) to illustrate their application and ver- 
satility. Subsequently, these theories are employed to predict the complete thermoelas- 
tic and strength behavior of boron-fiber - aluminum-matrix unidirectional composites. 
Some typical thermoelastic properties of off-axis and angle-ply composites are shown as 
a guide to design. Failure envelopes under off-axis loadings and under combined stress 
are presented. The effects of varying the transverse tensile strength and the intra- 
laminar shear strength on off-axis failure envelopes are investigated. The methods de- 
scribed in reference 21 are employed to predict some of the other important design fac- 
tors mentioned previously (such as transverse shear delamination criterion and buckling 
coefficient). The results are presented in graphical form. 

The final form of the equations used to generate the results is given in the appen- 
dixes. The physical considerations, assumptions, and detailed derivations are omitted 
here since they are available in the references cited. The experimental results of ref- 
ence 22 were used to evaluate some of the correlation coefficients in the semiempirical 
theories. 

The majority of the results presented herein are theoretical, and they were obtained 
assuming linear or approximately linear UFC behavior to failure. These results show 
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the application of composite mechanics in the characterization and generation of design 
data for fiber-reinforced metals. In addition, they bring into focus the type of informa- 
tion needed in design, the areas of difficulty in obtaining this information, the experi- 
mental data which are still needed, and the possible areas of immediate and future re- 
search. The discussion is supplemented by pertinent references in all cases. To this 
end, the detailed discussion of one specific composite system is sufficient. 

SYMBOLS 

A, A constants, eqs. (A 14) and (B5) 

a, a^ plate dimension; constants, eq. (B8) 

B constant, eq. (A14) 

C f, C m constants, eq. (A15) 

Cf, C constants, eq. (A16) 

Cj, CJn constants, eq. (A20) 

Cf*’ Cm constants, eq. (A21) 

d diameter 

E normal modulus, or elastic coefficients, or matrix elements 

F failure criterion function 

f functional 

G shear modulus 

H c heat capacity 

K, K bulk modulus or heat conductivity 

12’ 12 elastic and correlation coefficients, respectively, eq. (B12) 

k, k apparent and actual volume ratio, respectively 

Nj, number of fibers per end and number of layers, respectively 

P cr buckling load 

p parameter, eq. (B4) 

R transformation matrix, eq. (B18) 

S simple strength, failure or limit stress 

AT temperature change 
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t 


thickness 


W 

x, y, z 
1, 2, 3 
a 


^kv» ^kl’ %2’ %3 

^fT, fC, mT, raC, 22T, 22C, 12S 

6 

e 

C 

e 


V 

p 

o 


{ } 

[] 

[] T 

tr 1 

Subscripts: 


weight 

structural (load) axes and coordinates thereof 
material axes (1 axis coincides with fiber direction) 
thermal coefficient of expansion 

correlation coefficients for ply thermoelastic properties, 
eqs. (A15), (A16), (A20), and (A21), respectively 

correlation coefficients for heat conductivities, eqs. 

(A27) and (A28) 

correlation coefficients for strength, eqs. (B7) to (Bll), 
respectively 

in situ end, or strand, or yarn, packing array param- 
eter, eq. (A4), fig. 18 

interfiber spacing (not center to center), eq. (A3), 
fig. 19 

strain 

buckling parameter, fig. 22 

angle between load and filament direction 

buckling and delamination parameters, figs. 23 and 24, 
respectively 

Poisson ratio 

weight density 

stress 

strain magnification factor 

1 column array 

array 

array transpose 
array inverse 


C 


interface bonding 
compression 
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c composite property 

D debonding 

f filament property 

i ply index 

l ply property 

m matrix property 

p limiting property 

R residual stress 

S shear 

T tension 

v void 

x, y, z structural (load) axes direction 

1, 2, 3 material axes directions (1 axis coincides with the fiber direction) 
a, i 3 T or C, tension or compression 

BRIEF REVIEW OF PHILOSOPHY OF SEMIEMPIRICAL THEORIES 

The fundamental concept underlying semiempirical theories is summarized in the 
next three sentences. A simply mathematical model is selected which is mechanisti- 
cally representative of the physical state of affairs. The model is rendered predictively 
correct by strategically incorporated theory-experiment correlation factors. The cor- 
relation coefficients are then selected from a composite with some fiber content and 
made from a particular fabrication process. Subsequently, the theory is used to pre- 
dict properties of composites of other fiber contents but the same fabrication process. 
The following example illustrates the concept. 

The failure stress of a unidirectional fiber composite takes the following functional 
form: 


S l “ f [( k > d ’ N > A )f, V ’ k m’ ( E ’ v ’ G ’ S ’ € p) f m ’ S B’ a R’ °] 


( 1 ) 


where S denotes failure or limit stress, k, d, N, and A represent volume ratio, fila- 
ment diameter, number of fibers per end, and array packing, respectively, E, u, G, 
and are the normal modulus, Poisson ratio, shear modulus, and limit strain, re- 
spectively; a denotes stress state due to applied load; the subscripts l , f, v, m, B, 
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and R denote ply, fiber, void, matrix, bond strength, and residual stress, respec- 
tively. Conversion of equation (1) into an explicit form, if possible, would require com- 
plex statistical mathematical formalisms. On the other hand, the ply longitudinal ten- 
sile failure stress (fig. 1) can be expressed simply by using a rule-of-mixtures 

model (eq. (B7), appendix B): 

S l 11T = S fT^fT k f + ^mT k m 

where the barred quantities indicate actual volume content, and and /3 m rp are the 

experiment-theory correlation factors which convert equation (1) to (2). The param- 
eters /3^rp and embody all the variables which appear in equation (1) but are ab- 

sent from equation (2). Another way to think of these parameters is as fabrication proc- 
ess variables. Variables such k y , A, Sg, and Cg are associated with a particular 
fabrication process. As a result, the correlation factors have to be evaluated for some 
particular composite system for each fabrication process (refs. 17 and 18). Subsequent 
application of the theories will illustrate evaluation of the theory-experiment correlation 
factors. 

The semiempirical theory described herein will handle composites exhibiting linear 
or approximate linear behavior. Composites containing ductile fibers exhibit some non- 
linear stress- strain behavior. The strength properties of these composites can be pre- 
dicted by the semiempirical theory in conjunction with the well known incremental load 
concepts. 


J mll 


E 


fll 


(2) 


APPLICATION OF SEMIEMPIRICAL THEORIES TO SOME EXISTING DATA 

The experimental data in figure 2 are those reported in reference 19 for the ply 
transverse modulus g 2 - The dashed line in this figure was obtained by applying the 
finite element method to a square array of filaments (ref. 19). The solid line is the 
result predicted by the semiempirical theory with values for the correlation coefficients 
and /3 as noted in the figure. The agreement between semiempirical theory and 
the experimental data is quite good for a large range of k^ values. If the residual 
stress and bond strength vary appreciably over some range of k^ values, then different 
values of the correlation coefficients might be needed in different portions of this range. 

The experimental data in figure 3 are those reported in reference 19 for ply trans- 
verse failure stress S^22 t* The dashed line in this figure was obtained using the finite 
element method and a square array assuming complete interfacial bond (ref. 19). The 
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Figure 2. - Comparisons for transverse modulus of boron-aluminum unidirectional fiber 
composites. Filament modulus, Ep 55xl0 6 psi; matrix modulus, E^ = 10xl0 6 psi; 
filament correlation coefficient, pp 1.0; matrix correlation coefficient, p m = (1.0/k m r 6 ; 
void volume ratio, k y = 0. 


O Experimental data (fig. 4 of ref. 19) 
Finite element theory, square array 



Figure 3. - Comparisons for transverse strength of boron -aluminum 
unidirectional fiber composites. Matrix tensile strength, S mT = 

18 000 psi (ref. 19); filament modulus, Ep 55xl0 6 psi (ref. 19); 
matrix modulus, E m = lOxlO 6 psi (ref. 19); filament correlation coef- 
ficient, pp 1.0; matrix correlation coefficient, p m = (1. 0/k m ) 1/6 ; 
void volume ratio, k y = 0. 


solid lines are those predicted by the semiempirical theory with correlation coefficient 
/?22 t ( ec b (®9)) values as noted in the figure. Two values are required for ^22T in 
this case: &22T = 1- 08 for the range 0 < k f « 0. 2 and /^t = 1* 0 for 0. 2 < k f 0. 5. 
Distinct values of the correlation coefficient in a certain k f range are more typical; 
this is also one of the basic hypotheses for using semiempirical theories. 

It might appear, at first glance from the results in figure 3, that the finite element 
method is far superior to the semiempirical method in this case. However, the finite 
element method is far too complex and time consuming as compared with equation (B9). 
The complexity of the finite element method increases exponentially when it is called 
on to handle an anisotropic in situ matrix, anisotropic fibers, and void effects. 

The experimental data in figure 4 are the room temperature data reported in refer- 
ence 20 for ply longitudinal tensile stress j^rp. The values of the correlation coef- 
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o 

5 ? 


Experimental data (ref. 20), table I ; 

room temperature only) 
Semiempirical theory 
Used to evaluate correlation 
coefficient 



Figure 4. - Comparisons for longitudinal tensile strength of boron- 
aluminum unidirectional fiber composites. Filament tensile 
strength, Sfj= 510 ksi; filament modulus, Ef = 59xl0 6 psi; matrix 
modulus, E m = 10x10° psi; void volume ratio, k v = 0. 


ficients /3^ T and /3 m rp entering equation (B7) and the experimental points selected to 
evaluate the correlation coefficients are noted in the figure. Here, it is seen that values 
of the correlation coefficients can be selected for good theory- experiment agreement 
through any desired points. 


THERMOELASTIC PROPERTIES OF BORON-ALUMINUM UFC 

The constituent properties employed to generate ply thermoelastic properties for 
composite characterization and design are given in table I. The experimental data used 
to select the correlation coefficients are shown in table II. The corresponding values 
of the correlation coefficients entering the various equations in appendix A are listed in 
table III. 


Thermal Properties 

Ply thermal properties (three conductivities, three coefficients of thermal expan- 
sion, and the heat capacity) as functions of fiber content and assuming zero voids are 
presented in figure 5. The results in this figure were obtained using the appropriate 
equations in appendix A in conjunction with the constituent properties of table I and the 
correlation coefficient values of table III. As can be seen in figure 5, the heat capacity 
increases as the fiber content increases. The remaining thermal properties decrease 
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TABLE I. - MEASURED VALUES OF 


CONSTITUENT PROPERTIES 


Properties 

Boron 

fiber 

6066 -T6 alu- 
minum matrix 

E, 10 6 psi 

60.0 

10 

v, ratio 

0.2 

0. 3 

a, 10”^ in. /in. /°F 

2.8 

12. 9 

K, Btu/hr-ft 2 -°F/in. 

22.3 

84 

H, Btu/lb/°F 

0. 31 

0.23 

p, lb/in. 2 

0.085 

0.98 

S, ksi 

460 

a 52 

€ mpC’ in '/ in - 


a 0. 0052 

€ mpS’ in -/ in ' 


a 0. 00905 


a Values at yield. 


TABLE n. - MEASURED VALUES 
OF COMPOSITE PROPERTIES 
(REF. 22) 

[Fiber volume ratio, 0. 5. ] 


Property 

Value, 

psi 

E m 

34. 2X10 6 

*7 22 

23.4X10 6 

S l 11T 

164X10 3 

s nic 

183 

22T 

13.9 

S J22C 

15.3 

S il2S 

15.6 
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Heat conductivities, Kqj 2 2 33* 


TABLE m. - SOME BASIC PARAMETERS 
AND CORRELATION COEFFICIENTS FOR 
THERMOELASTIC PROPERTIES 


Property 

Value 

Basic parameters 

N f 

1 

d f 

0.004 


0 

k f 

0.35 < k f < 0.75 

Correlation coefficients 

% 

1.0 


4.0 


a 1.0 


a 2. 0 


a 1.0 


a 2. 0 

£ f 

a 1.0 


a 2. 0 

^kv 

a 1.0 

%i 

a 1.0 

%2 

a 1.0 

%3 

a 1.0 


a Estimate. 


80 i— 


60 h 


40 h 


20 




Figures. - Boron-aluminum unidirectional thermal properties. Void 
volume ratio, k v = 0. 
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with increasing fiber content. The literature is virtually void of measured data of 
thermal properties of fiber-reinforced metals. 


Elastic Properties 


Ply elastic properties (normal moduli, shear moduli, and Poisson ratios) are 
plotted in figure 6 as functions of fiber volume ratio and assuming zero void content. 


\n 

Q. 


O 

r*S 

cS* 




Figure 6. - Bororvaluminum unidirectional elastic properties. Void 
volume ratio, k v = 0. 


The results in this figure were obtained from the appropriate equations in appendix A 
using the constituent properties given in table I, the measured composite properties 
presented in table II, and the correlation coefficients of table III. It can be seen in fig- 
ure 6 that the moduli are increasing functions of fiber content while the Poisson ratios 
are decreasing functions. Note, particularly, that the transverse moduli E f 33 
increase rapidly with increasing fiber content. 

An extensive amount of experimental data is available for the longitudinal modulus 
n , and these data are in agreement with rule of mixtures predictions. Measured 
data for the remaining elastic constants are meager and more organized measured data 
are needed. 
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PLY STRENGTHS 


Ply strength, in the sense used here, denotes that stress level in the ply which 
causes incipient material nonlinear response (ref. 18). 

The ply strengths are characterized by the ply uniaxial strengths and a combined 
stress strength criterion (ref. 18). The constituent materials, the properties given in 
table I, the measured data in table n, and the correlation coefficients in table IV were 
used in conjunction with the equations in appendix B to generate ply strengths. 


TABLE IV. - CORRELATION COEFFICIENTS 
FOR STRENGTH PROPERTIES 
[Application limit, 0. 35 < < 0. 75. ] 


Correlation 

coefficient 

Value 

Correlation 

coefficient 

Value 

Pf T 

0.58 

a i 

8. 33 

%C 

1.00 

a 2 

52 000 

/3 mT 

^mC 

o o 
o o 

K 'l 12TT 
K 'l 12CT 

a 1.0 

a 1.0 

^22T 

0.313 

K 1 12TC 

a 1.0 

ft 22C 

0.343 

K l 12CC 

a 1.0 

&12S 

0. 462 




a Estimate. 


Uniaxial Strengths 

Five ply strengths (S mT , S znc , S z22t , S z22c , and S zl2g ) are required to char- 
acterize the uniaxial strength behavior of a ply (fig. 1). The graphs of these strengths 
as functions of fiber content when assuming zero voids are shown in figure 7. Measured 
values from table II are also shown in this figure. The curve labeled S z in figure 7 

limits the longitudinal compressive stress according to the debonding model (ref. 18). 

The results in figure 7 show that the S Z11C is matrix controlled for k f < 0. 5 and that it 
is debonding controlled for k f > 0. 5; S Z22T , S z22c , and S zl2g are relatively insensi- 
tive to fiber content in the range investigated. This is in contrast to the behavior of 
these same limit stresses which are quite sensitive to fiber content for nonmetallic ma- 
trix composites. 

Ply strengths for a ply with 5 percent void volume content are plotted in figure 8 
using the correlation coefficients from table IV since the theory accounts directly for 
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Figure 8. - Boron-aluminum uniaxial failure strengths. Void volume ratio, k v = 0. 5. 


the void effects. The results in figure 8 show that S^ 2 2 t> s i22C’ and S Z12S 

are slightly sensitive to fiber content at high fiber packing densities. The results also 
show that the ply longitudinal compressive strength is limited by debonding in plies with 
5 percent void content. 

Comparing similar curves in figures 7 and 8, it is readily seen that the presence of 


voids (porosity) has deleterious effects on ply strengths - particulary 
S 122 C’ and S U2S- 


S Z22T’ 


Researchers in the field of fiber- reinforced metals have investigated extensively the 


uniaxial longitudinal tensile strength for various fiber metallic matrix systems. How- 


ever, measured values of the other uniaxial strength properties are either rare or non- 
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existent. Measured data for the uniaxial strengths S^22T’ ^Z22C’ anc * ^Z12S 

are mandatory if fiber- reinforced metals are to become serious contenders for structural 
components and if their interfacial characteristics are to be better understood. 


Ply Combined-Stress Strength Criterion 

The combined- stress strength criterion for the ply used herein is described in 
equation (B12). The coupling coefficient ^ defined by equation (B13) and values 
for the correlation coefficients ^ a ’ @ = ^ are Si yen in table IV. 



Figure 9. - Boron aluminum combined-stress strength coefficient 
(applicable for 0 < k y < 0. 20). 


A plot of the coupling coefficient as a function of fiber volume ratio is given 

in figure 9. This coefficient is relatively insensitive to void content. The ply combined- 
stress strength criterion for a ply of 50 percent fiber volume content and zero voids is 
plotted in figure 10. Note that the 0^22 stress scale is different from the cx^ ^ scale 
in this figure. This causes the graph generated by the criterion to be distorted although 
it can be read more accurately this way. The various contours in the figure represent 
different levels of shear stress as noted. Measured values of the uniaxial ply strengths 
from table II are also shown in figure 10. 

The results in figure 10 show that the ply combined- stress strength is not sensitive 
to small values of the ratio a ^ This means that in testing for uniaxial normal 

limit strengths, a small amount of shear stress has negligible effects on the results. 

Failure envelopes for other fiber contents can be constructed by substituting corre- 
sponding ply uniaxial strengths from figure 7 and coupling coefficient values from fig- 
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Figure 10. - Failure envelopes under combined -stress for boron - 
aluminum. Fiber content, kf = 0.5; void volume ratio, k v = 0. 


ure 9 in equation (B12). Figure 8 can be used, or an analogous one, for eases in which 
the plies contain voids. 

There are no experimental data available for plies or unidirectional composites sub- 
jected to combined stress. For this reason, the correlation coefficient * s as " 

sumed to be one. Strategic experimental work in this area is needed for two important 
reasons: (1) to verify the failure criterion of equation (B12) or to provide insight and 
a basis for its modification, and (2) to generate measured data from which the correla- 
tion coefficient i2af3 can eva l ua -t e d- 


OFF-AXIS COMPOSITES 

Off-axis composites are unidirectional composites whose fiber direction is offset 
from the structural axis by a specified angle. The study of thermoelastic and strength 
behavior of off-axis composites serves two purposes: (1) generation of design data and 
(2) providing of additional insight into the composite structural behavior. 


Plies or UFC Subjected to Off-Axis Loadings 

Off-axis thermoelastic and strength characterization of a composite is readily de- 
termined from its characterization along the ply’s orthotropic axes using suitable 
transformations such as equations (B16) to (B20) (see appendix B). 
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Off-Axis Ply Thermoelastic Properties 

The strain- stress-thermal relations for an off-axis ply (assumed a thin plate) are 
given by 


or 


r ^ 

e cxx 


^ c cyyf 


\! cxy s 


1 

y cyx 

y csx 

E 

cxx 

E cyy 

G cxy 

^cxy 

1 

_ y csy 

E 

cxx 

E 

cyy 

G cxy 

v cxs 

_ ^cys 

1 

E 

L cxx 

E 

cyy 

G cxy. 


~\ 

cxx 


<p 


cyy 


> + at< 


txx 


cr 


a 


cyy 


> 


a 




cxy ) 


f € CX> = + AT < a c.yl 


cx J 1 cx 


cx j 


(3) 


where e, a, and AT are the strain, stress, and the temperature above some reference 

value; and E, G, v, and a are the normal modulus, shear modulus, Poisson ratio, 

and thermal coefficient of expansion, respectively. Note the Poisson ratios v„ cv and 

c sx 

v . These ratios couple in-plane normal and shear responses. Note also the thermal 
c sy 

coefficient of expansion a» „. Coefficients such as y , etc., and a; are present 

cxy csx cxy 

only in materials exhibiting anisotropic behavior. In equation (3), the following relations 
hold: 


y cyx _ y cxy . y cxs _ y csx . 

EE EG 

cyy cxx cxx cxy 

The matrices [E ] and { cn cx } in equation (3) are derived from the transformation 

equations (B 16) to (B20). Plots for the coefficients o'—, a , and a» , as functions 

cxx cyy cxy 

of the angle and for a particular fiber content (kj = 0. 5, 1^ = 0), are shown in figure 11. 

Analogous plots for the coefficients E, G, and v are shown in figure 12. The geometry 

of the off-axis ply and the ranges of the angle are also shown in these figures. The 

points to be noted in figures 11 and 12 are as follows: 

(1) The anisotropic coefficients a , u , and u vanish when the ply material 

cxy csx cby 

axes (axes 1-2) coincide with the ply structural axes (axes x-y). 


y cys _ y csy 

E G 
cyy cxy 
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(2) The anisotropic coefficients except the shear modulus attain significant magni- 
tudes at a certain angle. 

At these angles, residual stresses and severe end effects can be induced. Residual 
stresses and end effects are the primary obstacles in theory- experiment correlation and 
in interpreting experimental data. 

Plots for other fiber contents can be obtained by using figures 5 and 6 in conjunction 
with equations (B16) to (B20). In the case of voids, plots analogous to those of figures 
5 and 6 need to be generated first from the equations in appendix A for the specific void 
content. Heat conductivities for off-axis plies can also be obtained by transformation. 
The ply heat capacity remains invariant. 

Again, strategic experimental data to check the transformation equations are 
needed. 


Off-Axis Ply Failure Envelopes 

Plies or UFC subjected to off-axis loadings possess four distinct failure envelopes - 
namely, normal tensile load, normal compressive load, positive shear load (tends to 
extend the fibers), and negative shear load (tends to compress the fiber) (ref. 18). 

These envelopes are generated using values from figures 7 and 9 and table IV in equa- 
tion (B14) for normal loadings and equation (B15) for shear loadings. 

Plots of equation (B14) as a function of the angle and for a specific fiber volume 
content (k^ = 0. 5, k^ = 0) are presented in figure 13. The filament and loading direc- 



Figure 13. - Boron-aluminum off-axis normal load failure envelopes. 
Actual fiber content, k f = 0.5; void volume ratio, k v = 0. 
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tions and the experimental strength points ( 9 = 0° and 90°) from table II are also shown 
in the figure. The results in figure 13 indicate a rapid drop in the ply strength as load 
direction deviates from the filament direction for both normal tensile and compressive 
loads. This clearly illustrates the importance of fiber alinement when testing or de- 
signing for longitudinal strength. 

Analogous plots for shear load (eq. (B15)) are presented in figure 14. The filament 
and loading directions and the experimental strength points from table IE are shown in 
the figure. The upper portion of the figure is for positive shear load (as noted) while 



Figure 14. - Bororhaluminum off-axis shear load failure envelopes. 
Actual fiber content, kf = 0. 5; void volume ratio, k v = 0. 


the lower part is for negative load. The results in figure 14 indicate that the failure 
envelopes for off-axis plies are insensitive to shear loads. Experimental verification 
of the last comment is still needed. 

A last comment to be made for ply off-axis failure envelopes is that for each point 
on the uniaxial strength curves in figures 13 and 14 there corresponds a unique point in 
figure 10. For actual design purposes the ply off-axis failure envelopes are redundant 
and rarely used. However, these envelopes provide an insight into composite strength 
behavior. This will become clear from the discussion in the next section. 
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Effects of Intralaminar Shear and Transverse Strengths on 
Ply Off-Axis Failure Envelopes 


The ply off-axis failure envelopes can be utilized to yield additional useful informa- 
tion about composite strength behavior. This is illustrated in figure 15 where the ply 
off-axis failure strength for normal load and k^ = 0. 5 is plotted as a function of the 
angle for various combinations of the intralaminar shear and transverse uniaxial 
strengths. The standard curve in figure 15 is the curve labeled 1. The corresponding 
measured values from table II are shown in the figure at 0=0° and 90°. The curves 
corresponding to intralaminar shear variations are labeled 2 , 3, 4, and 5. The curves 
corresponding to transverse strength variations are labeled 6, 7, and 8. 

Three important observations can be made using figure 15: 

(1) The effect of fiber misalinement can be minimized with large increases in ply 
intralaminar shear strength. This finding is encouraging with respect to short fiber 
and whisker reinforced composites. 


Experimental data (ref. 22) 
Curve S l22T , S U2S , 

ksi ksi 

1 14 15. 6 

2 14 20. 0 

3 14 30.0 

4 14 40. 6 

5 14 60. 0 

6 20 15. 6 

7 40 15. 6 

8 60 15. 6 



Figure 15. - Off-axis norma! load failure envelopes for various values 
of s 122T- and %2S boron-aluminum composites. Actual fiber 
content, k f = 0.5; void volume ratio, k y = 0. 
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(2) The ply off-axis failure envelope is sensitive to transverse strength in the 
30° < 9 < 90° range. 

(3) The ply in angle-ply composites has four predominant failure modes: these are 
longitudinal tension or compression, 0 ° < 0 < 3°; intralaminar shear, 3° < 8 < 30°; 
shear or transverse tension or compression, 30° < 6 < 60°; transverse tension or com- 
pression, 60° < 8 < 90°. These observations stem from the fact that a ply within an 
angle-ply composite is reinforced by its adjacent plies, relative to intralaminar shear 
and transverse failure stress modes. This type of reinforcement has the same effects on 
the ply off-axis failure envelope as do the variations of intralaminar shear and transverse 
strengths noted in figure 15. 

Appropriate experimental data are lacking for ply failure envelopes for fiber- 
reinforced metals of potential use for structural components. Also, experimental data 
to verify the three points previously mentioned are needed. 


SOME PROPERTIES OF ANGLE-PLY COMPOSITES 

Properties of angle-ply composites are determined theoretically by laminate 
theories (refs. 24 to 26). In general, lamination theory yields the strain-stress thermal 
relations, axial stiffnesses, coupling stiffnesses, bending stiffnesses, and thermal 
forces. Herein, only results for the thermal coefficients of expansion and the elastic 
constants of a symmetric (in-plane shear and bending) angle-ply composite are pre- 
sented. The results are for a composite with 50 percent fiber content, zero void con- 
tent, and the stacking sequence (6, -8, - 8 , 8). The results were obtained using the con- 
stituent properties in table I and equations (B2 1) and (B22) . 

Plots of the thermal coefficients of expansion ce and a as a function of the 

cxx cyy 

ply angle are shown in figure 16. Analogous plots for the elastic constants are shown 



Figure 16. -Thermal coefficients of expansion fair boron-aluminum 
angle-ply composites. Actual fiber content, = 0. 5; void volume 
ratio, k v = 0. 
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Figure 17. - Elastic constants for_boron-aluminum angle-ply com- 
posites. Actual fiber content, k f * 0. 5; void volume ratio, k v - 0. 


in figure 17. Angle-ply composites with symmetric stacking sequence relative to in- 
plane shear do not have coefficients such as a^y and v csx , etc. The results in fig- 
ure 17 indicate that the shear modulus of symmetric angle-ply composites, with low 

(E /E ) ratio, is relatively insensitive to ply orientation angle, 
cxx cyy 

Particular care should be taken in fabricating angle-ply composites. Factors such 
as ply stacking sequence, slight fiber misalinements in the ply, slight ply misorienta- 
tion from the structural axes, and variation in fiber and void content from ply to ply 
result in several distinct coupling responses of the type described in reference 27. 
Here, suffice it to say, that all these coupling responses affect the structural response 
of angle-ply composites. This makes it difficult (if not impossible) to interpret exper- 
imental results and to correlate theories with experimental data. 

Experimental data for angle-ply fiber-reinforced metallic -composites is limited. 
Strategic experiments to obtain some measure of the various coupling response effects 
on composite structural behavior will prove very beneficial. 


ADDITIONAL ASPECTS OF COMPOSITE CHARACTERIZATION MECHANICS 

Several other aspects of composite characterization mechanics, which have not 
been dealt with specifically herein, merit mentioning. Some others need verification, 
modification, or extension, and then there are those yet to be developed. Some of these 
are mechanics for constituent load transfer at the interface, composite toughness, com- 
posite internal damping, lamination edge and internal residual stresses, geometry and 
material nonlinearities (including plastic behavior), temperature and time effects in- 
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eluding creep, and cyclic loading and fatigue. References 28 to 37 should provide some 
direction in all these areas. 

A deeper assessment and carefully formulated programs (both analytical and ex- 
perimental) to investigate various aspects of the previously mentioned areas will be of 
great service to the fiber-reinforced metals research and design community. 


OTHER DESIGN FACTORS 

Composite thickness, interfiber spacing, composite density, delamination, and 
buckling parameters are important variables in component design. Some of these vari- 
ables can be related to fiber geometry, packing array, void content, and other constit- 
uent properties by elementary relations described in reference 21. Results are pre- 
sented here for a boron aluminum composite. 

The ply thickness is plotted in figure 18 as a function of the actual fiber content and 



Figure 18. - Boron-aluminum ply thickness. Packing array parameter, 
= 1.0; void volume ratio, k v * 0. 


for a square array. The curve in this figure was generated from equation (A 4) with 
Nj = 1, /3 t = 1 (the parameter can be selected from measured data for different ar- 
rays; ref. 20, and d^ = 0. 004 inch. Analogous plots for boron fibers with different 
diameters can readily be generated from equation (A4) The curve in figure 18 indicates 
that the ply thickness decreases rapidly with increasing fiber content. This is significant 
because the buckling resistance of a structural component varies with the thickness cubed. 
Note the presence of voids increases the ply thickness, if the apparent fiber volume re- 
mains constant. 

A plot of the interfiber spacing against actual fiber content and for a square array 
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Figure 19. - Boron-aluminum interfiber spacing. Void volume ratio, 
k v =0. 


is shown in figure 19. The curve in this figure was generated from equation (A3) with 
= 0. 004 inch. Analogous plots with different fiber diameters can be readily generated 
from equation (A3). The results in figure 19 show a rapid decrease in the interfiber 
spacing particularly at high packing densities. The designer of fiber composites should 
keep the curve in figure 19 in mind. This curve is important because transverse 
strength, intralaminar shear strength, and fiber-matrix load transfer depend strongly 
on the interfiber spacing. 

The composite weight density p c is plotted against actual fiber content k^ in fig- 
ure 20. As can be seen in the figure the composite density is a linearly decreasing func- 
tion of the actual fiber content for the boron- aluminum composite. The presence of voids 
in a composite reduces its weight density if the actual fiber content is held constant. 



Figure 20. - Boron -aluminum composite density. Void 
volume ratio, k v = 0. 
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Figure 21. - Boron-aluminum ply weight. Void volume 
ratio, k v = 0. 
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Figure 22. - Boron-aluminum buckling parameter for 
simply supported square plates. Void volume ratio, 
k v =0. 


An interesting parameter in trade-off and efficiency studies is the ply weight per 
unit area. This parameter is plotted against k f in figure 21 for square array packing. 
The resulting curve for this case shows that the ply weight decreases approximately in 
a hyperbolic fashion with increasing k^. 

The buckling resistance of a structural component is another parameter commonly 
employed in trade-off and efficiency studies. The bucling parameter £ for a square 
simply supported plate is plotted against k^ in figure 22. The plate in this case is 
loaded along the fiber direction and is made from unidirectional material with square 
array packing and a fixed number of layers (ref. 21). The plot in this figure shows a 
very rapid decrease of the buckling load at low fiber contents and a mild decrease at 
k^ > 0. 3 with increasing fiber content. The decrease of the buckling parameter at the 
low fiber volume contents comes from the fact that buckling resistance varies with the 
thickness cubed while it varies linearly with modulus. As is illustrated in figure 18, the 
thickness decreases as the fiber content increases for the square array fiber packing. 
The important point to be kept in mind from this discussion is that low fiber content 
composites are more efficient in resisting buckling when the number of layers is the 
limiting condition. Analogous results when the weight is held constant are plotted in 
figure 23. The curve in this figure shows that the buckling load is approximately a 
parabolic increasing function of k f and does not possess an optimum. The conclusion 
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Figure 23. - Boron-aluminum weight buckling parameter for 
simply supported square plate. Void volume ratio, k v = 0. 



Actual fiber content, k f 


Figure 24. - Boron-aluminum transverse shear delamination parameter. 
Void volume ratio, k v = 0. 
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in this case is that high fiber content composites resist buckling loads more efficiently 
when the composite weight is the limiting condition. 

Ply delamination, either due to adjacent ply relative rotation (scissoring effect) or 
interply transverse shear (through the thickness, horizontal), is a possible failure mode 
for multilayered composites (ref. 21). The delamination parameter for transverse 
shear is plotted against k^ in figure 24. The equation of the curve and the correspond- 
ing geometry are also shown in the figure. The notation in the equation is as follows: 

e„„_, interply transverse shear strain; X delamination parameter; cr _, interply 
cxz s cxz 

transverse shear stress; G^ xz , interply shear modulus in the x-z plane (ref. 21). Re- 
place x with y for the equation in the y-z plane. The resulting curve in figure 24 indi- 
cates that the delamination parameter for transverse shear is very sensitive to k^ with 
a maximum k^ — 0.65. This means that low and very high fiber content boron- aluminum 
composites are less susceptible to delamination than are those of intermediate k^ values. 
Results for the delamination parameter due to adjacent ply relative rotation are not pre- 
sented because evaluation of this parameter requires specific measured data which are 
not available at this time. Reference 21 describes the type of data needed. Measured 
data for several of these factors and theory- experiment correlation are still needed. 


CONCLUDING REMARKS 

The semiempirical theories developed primarily for nonmetallic composites can 
also be employed in characterization and design mechanics for metallic composites 
which exhibit linear or approximately linear behavior. This is important in optimum 
designs of components made from these materials where the fiber content is taken as 
a design variable. It is also important in designing experiments and interpreting experi- 
mental results. 

The results of this investigation lead to the following conclusions concerning boron 
fiber - aluminum matrix composites: 

1. The uniaxial strengths S^ 22 t> ®I22C’ anc * 12S are relatively insen- 

sitive to fiber content in the range investigated when voids are absent but are sensitive 
at high fiber contents in the presence of voids. 

2. The longitudinal compressive strength is limited by both matrix yield (k^ < 0. 5) 
and debonding (kj > 0. 5) in the absence of voids. This strength is governed by debonding 
(0. 40 < kf < 0. 75) when the void content is about 5 percent. 

3. The failure envelopes for off-axis normal loadings show rapid loss of strength 
when the angle between fiber and load directions increases in the 1° < 9 < 30° range. 
The strength remains practically invariant at angles greater than 30°. 

4. The failure envelopes for off-axis shear loadings are practically invariant with 
angle. 
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5. The normal tensile load failure envelope is quite sensitive to variations in intra- 
laminar shear strength for load angles less than 30° and to variations in transverse ten- 
sile strength for load angles greater than 45°. 

6. The predominant failure modes for off-axis normally loaded composites are 
longitudinal tension or compression when 0° < 9 -= 3°, intralaminar shear when 

3° < 0 <= 30°, intralaminar shear or transverse tension or compression when 
30° < 9 ?S 60°, and transverse tension or compression when 60° < 9 < 90°. 

7. The shear modulus of off-axis and angle-ply composites is practically invariant 
with angle at low E £n/ E £22 ratios * 

8. The ply thickness, interfiber spacing, composite density, and ply weight per unit 
area decrease with increasing fiber content. 

9. The buckling load of a simply supported square plate made from a unidirectional 
boron aluminum composite, in which fibers are arranged in a square array and which is 
loaded in the plane of the plate, decreases with increasing fiber content when the number 
of plies is held constant but increases when the weight is held constant. 

10. The delamination criterion indicates that the transverse shear strain in the ma- 
trix increases with fiber content up to about 65 percent fiber content and then it de- 
creases very rapidly. 

11. Experimental data for thermal properties, transverse strength properties, and 
combined failure criteria are needed. Well defined programs to investigate tempera- 
ture, time, material, and geometry nonlinearities, damping, impact, and cyclic load- 
ings are required. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, February 2, 1970, 

129-03. 
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APPENDIX A 


EQUATIONS FOR PLY GEOMETRIC, ELASTIC, AND THERMAL PROPERTIES 

The equations of the semiempirical theory are given herein. Most of the detailed 
derivations are rather lengthy and are not included here; however, they are given in 
appendixes A and B of reference 26. It is assumed that the various thermoelastic prop- 
erties of the constituents in their precomposite state are known. It is further assumed 
that the apparent filament volume ratio kj is given. When the axiom that the volume 
(weight) of the whole equals the sum of its parts is used, it can be shown that 


k f = (1 - Vk f 

(Al) 

k m = d - k f )(l - v 

(A2) 


where the barred quantities denote actual volume ratio. The interfiber spacing 6 and 
the ply thickness t^ are, respectively, given by 

6 = 


and 


\W 2 

where ^ is the in situ array rectangle length to width ratio (ref. 21). Let the normal 
strain- stress - temperature relations along the orthotropic axes for ply, filament, and 
matrix be given, respectively, by 



{e t } =[E l ]{a l } + AT{o- z } 

(A5) 

{e f } - [Ej]{ cr f } +AT(« { } 

(A6) 

{ e m^ = [Em-H cr m^ +AT{a' m } 

(A7) 
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The normal elastic properties and the thermal coefficients of expansion are given, re- 
spectively, by 


where 


[ E J - [CfJ T [ E f][Cf Z ]k f + [C mZ ] T [E m ][C m Jk m 

<«z> = [C fz ] T { o,f } k f + [C mi ] T { } k m 


^lfj 



*lli) 

1 


r *l 11) 

< °Z22 | 

S 


€ l 22 

>; 


■ 

a l 22 > 

23> 



^ e Z33> 



<- a l 33-> 



1 

_ v l 21 

_ v l 31 " 




E Z 11 

E Z 22 

E Z 33 



_ v l 12 

1 


_ v l 32 


E Z11 

E Z 22 

E Z 33 




_ v l 13 

v l 23 

1 




. E Z11 

E Z 22 

E Z 33 . 



(A8) 

(A9) 


(A10) 


(All) 


Analogous equations for filament and matrix are obtained by replacing the subscripts in 
equations (A10) and (All). 

The filament and matrix stresses ply stress relations are 


r„ 
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(A12) 
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The parameters ^ and /3 are chosen such that theory and experiment correlate as 
well as possible. When computing the ply thermal coefficients of expansion from equa- 
tion (A9) use 
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J 


instead of equation (A1 5). Here /3^ and /3 are chosen such that theory and experi- 
ment correlate for the thermal coefficients of expansion. 

The ply shear moduli are given by 


G Z12 


T ml2 


hi 12 

ri’2p 

W Cr fl2 


kj + 
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C ' 2 
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13 


G ml3 


m13 k f 


pi * 2 r 

u f ^f 13 



where 



(A17) 


(A 18) 


(A19) 


(A20) 


(A21) 


The parameters /3j, /3j’, and /3^ are empirical factors and are to be chosen such 

that predicted and experimental results for the shear moduli correlate. The empirical 
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factors /3jl and (3' m can be selected from horizontal-shear-beam tests where the fila- 
ments are parallel to the beam axis. The empirical factors j8£' and can be se- 
lected from horizontal- shear -beam tests where the filaments are normal to the beam 
axis. It is noted that the condition 


E £ 33 

2 ( 1 + V 132 ) 


22 

2 ^ 1 + 


(A22) 


should be satisfied when the ply possesses transverse isotropy in the 2-3 plane. For 
this case, 


E 


7 33 


= E 


l 22 


v l 32 “ V 12Z 


The hydrostatic ply bulk modulus is given by 





^12 : ^23 
E Z11 E 7 22 


+ ^3\ 

E l 11 / 


The ply heat capacity is given by 


% = [j-J [ H cf p A + Wm] 


(A23) 


(A24) 


(A25) 


where is the ply weight density given by 

p l = [ P A + p m^m] ^ A26 ^ 

The heat conductivity of the matrix containing randomly distributed spherical voids is 
given by 


K = K 
maa m ota 


2 ^kv^maa! + E v 




K v ) 




(A27) 
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where denotes the air conductivity and the parameter j3y. v is the empirical factor 
chosen so that predicted and experimental results correlate. The subscripts aa re- 
fer to the orthotropic axes of the matrix. The heat conductivities of the ply are given 
by 


K l 11 " %l k f K fll + k m K mll 





K l 22 “ K m22 

1 " %2 V 


f i I < 22 ) 



%2-^f \ K f 22 / _ 


(A28) 


33 “ K m33 





It is important to note that k m is used in the first of equations (A28) instead of k m 

since K includes the void effects. The parameters /3, 1; etc. are semiempirical 

m 01 ot iv x 

factors chosen so that predicted and experimental results correlate. 


35 



APPENDIX B 


EQUATIONS FOR FAILURE STRESSES (STRENGTHS) AND TRANSFORMATION 
Matrix-Strain Magnification Factors and Void Effects 

The transverse and shear matrix-strain magnification factors (appendix of ref. 18 
or ch. 2 of ref. 26) are 


<P 


p22 


-[—l 

L 1 + p( a 


p(A - 1)J 
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where 


P = 



and 


A = 1 " / E m22 \ 

1 ' l ml2 V m21\ E f22/ 

The maximum void effect factor is given by 


(B3) 


(B4) 


(B5) 
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(B6) 
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The parameters E, G, v, and a denote longitudinal modulus, shear modulus, Poisson 
ratio, and stress, respectively. The subscripts v, f, and m denote void, filament, 
and matrix, respectively, and the numerical subscripts correspond to the filament di- 
rections depicted in figure 1. The variable k y is the void content, and k f and k m 
are defined by equations (Al) and (A2), respectively. 


Uniaxial Strengths 

The uniaxial strengths 22T’ ^£22 C j and are given by 

(ref. 18) 
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where S fT , e mpT ? e mpC ? anc * e mpS re P resent filament bundle strength, matrix 

compressive strength, matrix tensile yield strain, matrix compressive yield strain, and 
matrix shear yield strain, respectively; E and G denote normal and shear moduli; k 
is the volume fraction, and <p and /3 y represent the strain magnification factors and 
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the maximum void effects, respectively; and the remaining /3's denote the simple 
strength correlation coefficients. 


Combined-Stress Strength Criterion 

The ply failure criterion for combined stress is given by (ref. 18) 


Y{a L , Sj , Kj i 2 ) - 1 - 
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7 22/3 
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1 12S 


n 12 s/ 
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(B12) 


where F denotes the combined- stress strength criterion as follows: 

F(ct^ , , Kj 12 ) >0 no failure 

F(a^ , ^ - 0 onset of failure 

F(a^ , ^ 0 failure condition exceeded 

and where subscripts a and /3 denote T (tension) or C (compression), the applied 
stress state determined from the stress analysis, S the unidirectional filamentary com- 
posite (UFC) simple strength either determined from the equations described previously 
or measured experimentally, and ^a/3 the theor y- experiment correlation coefficient 
determined as described in reference 18. The coefficient K ^ ^ 2 is given by 

K (! + 4 v ll2 - ^ 13 )E g2 2 + (*• ~ 23**7 11 (B13) 

i 12 r ”1 i/o 

[ E £ 11 E Z 22^ 2 + v l\2 + v l 13^ 2 + v l 21 + v l 23 J 

where and v l denote UFC modulus and Poisson's ratio, respectively. The sub- 
scripts 11, 13, etc. refer to the corresponding axes in figure 1. For the case of iso- 
tropic material, equation (B13) reduces to one, as can be verified by direct substitution, 
and equation (B12) reduces to the well known von Mises criterion. 

Off-axis failure e nvelopes fo r UFC. - There are four possible off-axis loadings on 
a ply - namely, normal tensile force, normal compressive force, positive shear force, 
and negative shear force (ref. 18). The corresponding equations are, respectively, 
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where is the failure strength, 0 the angle between filament and load directions, and 
a either T (tension-tensile load envelope) or C (compression-compressive load enve- 
lope). When a = T and j3 = C, equation (B15) describes the positive shear envelope; 
when a = C and 0 = T, equation (B15) describes the negative shear envelope. 


Transformation Equations for Off-Axis UFC 


The two-dimensional elastic and thermal constants of off-axis UFC are given by 
(refs. 23 and 24) 
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and where 9 is the angle between the x-axis and the fiber direction. The elements in 
[E cx ] anc * ^ a cx } are obtained by replacing l by cx in equations (B19) and (B20). 


Equations for Angle-Ply Composites 

Let z^ and z^ represent distances from some reference plane to the top and bot- 

i-T- 

tom of the l ply, respectively. Then, neglecting interply layer contributions, the 
elastic and thermal coefficients of expansion of angle-ply composites are given by, re- 
spectively, 
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i“cx) 'I E cx] L t R ll) T [ E ur l (aji) 
i=l 


(B22) 


' I 


(B23) 
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where [R^J, [E^j], and {or^} are matrices representing transformation (eq. (B18)), 
elastic constants (eq. (B19)), and thermal coefficients of expansion (eq. (B20)) for the 

|j. 

i ply, respectively. 
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